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We use ab initio density functional calculations to study hydrogen- induced disintegration of single- 
and multi-wall carbon fullerenes and nanotubes. Our results indicate that hydrogen atoms prefer- 
entially chemisorb along lines in sp^ bonded carbon nanostructures, locally weakening the carbon 
bonds and releasing stress. For particular structural arrangements, hydrogen helps to relieve the 
accumulated stress by inducing step-wise local cleavage leading to disintegration of the outermost 
wall. 
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Exposure to hydrogen is known to cause embrittle- 
ment in metals[l[, a major concern in materials sci- 
ence. It is an intriguing question, whether local deco- 
hesion due to chemisorbed hydrogen may cause similar 
damage in sp^ bonded carbon nanostructures, such as 
fullerenes and nanotubes [2|, and whether it may have 
helped to efficiently cleave graphene[2i]- So far, past re- 
search has focussed on beneficial side-effects of the in- 
teraction of hydrogen with nanocarbons, including eas- 
ier separation/debundling and the prospect of reversible 
hydrogen storage [3, Since both fullerenes and nan- 
otubes are under internal stress due to the local curva- 
ture, hydrogen-induced local decohesion may facilitate 
fracture and cause irreversible damage in nanocarbons 
as it does in metals [1]. 

Here we present an ab initio density functional study of 
bonding and structural changes in sp^ carbon nanostruc- 
tures exposed to hydrogen. Our results indicate that hy- 
drogen atoms preferentially chemisorb along lines on the 
outermost wall of single- and multi-wall carbon fullerenes 
and nanotubes, locally weakening the carbon bonds and 
releasing stress. For particular arrangements of adsorbed 
hydrogen atoms, we identify the microscopic process of 
step-wise local cleavage leading to the disintegration of 
the outermost wall as a way to release the accumulated 
stress. 

To gain fundamental insight into the effect of hydro- 
gen on the structural integrity of sp^ bonded carbon 
nanostructures, we studied the total energy and opti- 
mum structure of pristine and hydrogenated fullerenes 
and nanotubes. Our results are based on the density 
functional theory (DFT) within the local density approx- 
imation (LDA) and ffist-principles pseudopotentials. We 
used primarily the SIESTA [s] code with a double-^ ba- 
sis set with polarization orbitals, the Perdew-Zunger[6j 
parametrization of the exchange-correlation functional, 
and norm-conserving ab initio pseudopotentials p?] in 
their fully separable form [8]. The accuracy of our results 
was verified using the PWSCF plane- wave code with ul- 
trasoft pseudopotentials. We used 10~^ eV/A as a strict 



gradient convergence criterion when determining opti- 
mum geometries. Constrained structure optimization has 
been used to investigate energy barriers, transition states, 
and generally to identify the optimum reaction pathway. 

Published theoretical results include studies of the 
bonding character of isolated H atoms on graphite 
flakes, nanotubes and fullerenes [1, [13, and the pref- 
erential arrangement of hydrogen pairs on a graphene 
layer (12I, ll3|, Other investigations address the rel- 

ative stability of hydrogen pairs for different adsorp- 
tion geometries jlH and coverages [l6| of hydrogen on 
graphene. Theoretical studies of hydrogen-covered nan- 
otubes focussed on adsorbate-related electronic struc- 
ture changes [17] and the relative stability of different 
adsorption arrangements jioi 1 151 . llSl . llOl . l20|. More re- 
cent studies considered the possibility of axial cleavage 
in nanotubes in the presence of exohedrally adsorbed 
H atoms [lol, llSl, IZ^- No results have been reported so 
far about the origin of hydrogen-induced structural dam- 
age in fullerene-like structures and its relation to stress- 
induced cleavage of nanotubes. 

Same as in graphene, each carbon atom in fullerenes 
and nanotubes has three neighbors to form very stable 
sp^ bonds. Unlike in planar graphene, the most sta- 
ble sp^ carbon allotrope, stress is associated with the 
nonzero local curvature in fullerenes and nanotubes. The 
small fullerenes C20 and Ceo are spherical, with all atoms 
equivalent and subject to the same stress. With in- 
creasing size, however, free-standing fullerenes start re- 
sembling icosahedra, with twenty near-planar triangular 
facets spanning the surface. In these structures, stress 
associated with higher local curvature concentrates near 
the thirty edges connecting adjacent pentagons. In car- 
bon nanotubes, on the other hand, all atoms are equiv- 
alent and the stress associated with local curvature is 
distributed uniformly throughout the structure. 

Even though fullerenes with more than 100 atoms have 
not been observed, large icosahedral shells are believed to 
form the outermost shell of mult i- wall fullerenes, called 
bucky onions [22]. Focussing on the faceted Cgeo hollow 
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FIG. 1: (Color online) Equilibrium structure and bonding 
character in H-covered fullerenes and nanotubes. (a) Hydro- 
gen chemisorption on the Cgeo fullerene along strained lines 
with larger curvature, (b) Lines of chemisorbed hydrogen, 
shown in end-on and side view, yield the most stable adsorp- 
tion on the (6,6) carbon nanotube. (c) Relaxed atomic struc- 
ture near a line of chemisorbed hydrogen atoms, (d) Differ- 
ence charge density Ap indicating hydrogen-induced sp^sp^ 
local bonding change in graphitic nanocarbons. Contours of 
Ap in the plane indicated in (c) are superimposed with the 
atomic lattice. 



shell structure in Fig. [TJa), increased stress makes edge 
atoms more reactive and prone to bind hydrogen. Molec- 
ular hydrogen does not dissociate on planar graphene, 
but does so with an energy gain on fullerenes, with the 
dissociated hydrogen pair preferentially binding on top 
of adjacent carbon atoms. In the following, we will de- 
fine adsorption energies per hydrogen pair, related to the 
free H2 molecule. As expected, the energy gain depends 
on the local stress in the substrate, with values ranging 
from 1.20 eV at the twelve pentagonal corners to 1.35 eV 
at the edge sites and 0.34 — 0.67 eV on the facets of the 
C960 fullerene. 

Calculations at higher coverages suggest that after 
the first two atoms are adsorbed, additional hydrogens 
should preferentially adsorb along straight lines on the 
substrate [10]. We believe that formation of this ener- 
getically preferred adsorption pattern can be achieved 
by surface diffusion of hydrogen atoms at tempera- 
tures found during many hydrogenation reactions, since 
the activation barrier for sigmatropic rearr ang ement of 
chemisorbed hydrogen atoms is only ^1 eV[lO]. As sug- 
gested earlier and illustrated in Fig.[T](c), carbon atoms at 
icosahedral edges in large fullerene cages are more reac- 
tive than the facet sites. Since the dissociative chemisorp- 
tion energy of hydrogen along the edges exceeds that on 
fullerene facets by ^0.7 eV per H2 molecule, we expect 
hydrogen atoms to adsorb preferentially along lines con- 
necting pentagon corners. 

Even though there are no preferential adsorption sites 
on a carbon nanotube, atomic hydrogen is believed to 
adsorb preferentially along lines parallel to the tube 



axis [10!, 1^, m to release stress. Possible initial adsorp- 
tion geometries for hydrogen on the (6, 6) carbon nan- 
otube are shown in Fig.^b). In the left panel, we depict 
the optimum nanotube geometry following hydrogen ad- 
sorption along three zigzag lines parallel to the axis. The 
calculated binding energy of 0.54 eV per hydrogen pair 
reflects also the energy associated with the tube relax- 
ation to a triangular cross-section. Presence of hydro- 
gen apparently locally reduces the flexural rigidity of the 
nanotube wall, converting its circular cross-section into 
a triangle. 

Other favorable local adsorption patterns include exo- 
or endohedrally adsorbed hydrogens forming double-lines 
along the nanotube axis. The initial geometry of such 
a combined exo/endohedral adsorption arrangement is 
shown in the middle panel, and an exohedral double-line 
arrangement of hydrogens is shown in the right panel of 
Fig-IHb). As we discuss later on, hydrogen-induced local 
decohesion causes significant structural rearrangements 
in both cases. The energy gain associated with hydrogen 
chemisorption was found to be 0.87 eV per hydrogen pair 
in this case. 

To obtain fundamental insight into the origin of 
hydrogen-induced local decohesion that is general and 
independent of a particular graphitic structure, we stud- 
ied the effect of a line of chemisorbed hydrogen atoms 
on the structure and bonding in a graphene layer, as 
shown in Fig. [D^c). Hydrogen atoms bind preferentially 
on top of carbon atoms, causing a local buckling. A line 
of chemisorbed hydrogen atoms is expected to form a 
crease in the initially planar graphene [3]. In Fig.[T](c) we 
constrained the crease angle to equal that between adja- 
cent facets on the icosahedron in Fig.^a). To get insight 
into the nature of the adsorption bond, we inspected the 
charge density difference Ap = p(H/graphene) — p(H) — 
p(graphene), representing the charge redistribution in the 
system upon hydrogen adsorption. The spatial distribu- 
tion of Ap(r) is shown in Fig. [TJd) in the plane schemat- 
ically indicated in Fig.llfc). Hydrogen adsorption causes 
only local charge rearrangement leading to a charge ac- 
cumulation in the H-C a bonds, depleting the charge in 
the local C-C bond region, as expected in case of sp^ 
bonding. Both the charge rearrangement and the local 
buckling are strong indicators for a local sp^^sp^ transi- 
tion. As suggested above, hydrogen pairs bind preferen- 
tially to neighboring carbon atoms, thereby changing an 
initially strong sp^ bond to a weaker, strained sp^ bond. 
This local decohesion resembles to some degree the mi- 
croscopic mechanism of hydrogen embrittlement in some 
metals [1]. 

Besides decohesion, local change from sp^ to sp^ bond- 
ing is also associated with a net increase in the C-C bond 
length by typically 8%. In partially hydrogen- covered 
systems, carbon atoms connected to lines of hydrogen are 
under stress. The consequences can be best illustrated 
in the fullerene structure shown in Fig. [ija), where hy- 
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drogen only attaches to the edges. The hydrogen- free 
planar facets can be thought as incompressible, same as 
small graphene flakes. Since the corners of the triangu- 
lar facets determine the edge length, sp^ bonded carbon 
atoms along the edges are prevented from increasing their 
interatomic distance, subjecting the edges to stress. The 
energy associated with accumulated stress increases lin- 
early with the edge length L. For a critical value Lcrit^ 
this energy will exceed the activation barrier AE to break 
a C-C bond as a first step to initiating a cleavage that 
eventually causes the disintegration of the structure. 

As we show later on, we find that AE^l.7 eV. This 
amount of stress energy is accumulated in a hydro- 
genated fullerene edge of critical length Lcrit^lS nm. We 
thus expect hydrogen-induced cleavage to occur only in 
fullerenes with L > 15 nm, with a diameter exceeding 
^28 nm and containing more than 10^ C atoms. Such 
large fullerenes never occur as single-wall structures, but 
may form the outermost wall of a multi-wall bucky onion. 
In large bucky onions, the spherical walls inside exert ra- 
dial pressure on the outermost wall, thus enhancing the 
hydrogen-induced stress and reducing the critical diam- 
eter for disintegration of the outermost wall. We expect 
that multi-wall fullerenes exposed to hydrogen should 
peel layer-by-layer from the outside in a process appear- 
ing as exfoliation. 

In the following, we study the microscopic mechanism 
leading to cleavage at the minimum energy cost. To 
obtain a general understanding of the fracture process 
that does not depend on a particular system, we created 
generic creases by placing a periodic array of hydrogen 
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FIG. 2: Zipper mechanism of hydrogen-induced disintegra- 
tion of strained sp^ carbon nanostructures. (a) Periodically 
warped graphene as model strained system, (b) Optimum 
structure and energetics for breaking a bond pair, (c) Se- 
quence of bond cleavages resulting in a crack. 



lines along the armchair direction on alternating sides 
of an infinite graphene monolayer, as seen in Fig. [Sfa). 
The energetics and atomic arrangement during the ini- 
tial cleavage are depicted in Figs.[2](b) and[2](c). Cleavage 
at the crease is caused by locally weakening the affected 
bonds. We found that this is best achieved by adsorbing 
additional hydrogen atom pairs on both sides of the ini- 
tial hydrogen line, as shown in the left panel of Fig.[2jc). 
Such a local arrangement of adsorbed hydrogens is a fa- 
vorable prerequisite to cleave an edge segment. We found 
that the energy needed to break one bond is larger than 
the energy to break two adjacent bonds, shown in a dark 
color in the left panel of Fig. [21(c). Continuing this con- 
certed bond breaking mechanism should cleave the edge, 
causing the structure to disintegrate. 

The energetics of cleavage in Fig. Efb) shows the total 
energy change with respect to the initial structure as a 
function of the bond length d across the crack, which we 
use as the reaction coordinate. We consider the full range 
of d values occurring during the cleavage reaction and, for 
each value of we globally optimize the structure. Con- 
sequently, the calculated transition paths for the disin- 
tegration of hydrogenated fullerenes and nanotubes rep- 
resent the optimum cleavage path in unconstrained sys- 
tems. Following an initial energy investment of 1.7 eV, 
the bonds break abruptly, leading to the final structure 
depicted in the right panel of Fig. [Sfc). There is a net 
energy gain of ^2 eV associated with the bond breaking, 
caused by releasing the accumulated stress. The initially 
stressed sp^ bonded crease segment converts into a pair of 
more stable sp^ bonded, hydrogen terminated graphene 
edges. 

We next calculated the energy barriers associated with 
breaking specific bonds in search of the optimum path- 
way to complete the crease cleavage after it was initiated. 
The optimum sequence of bonds to break, along with the 
intermediate structures and the energetics of the process, 
is depicted in Fig.[2jd). Interestingly, the smallest energy 
investment to propagate the fracture following the first 
step in Fig. EJ^b) does not involve the neighboring pair of 
C-C bonds on the opposite side of the ridge, but rather 
the more distant pair of C-C bonds on the same side of 
the ridge. Our results indicate that the energy barrier 
to propagate the fracture decreases at successive steps. 
As a matter of fact, it may be reasonable to assume that 
the energy released following the initial step may acti- 
vate the zipper-like cleavage at the crease without fur- 
ther energy investment. Once cleavage is initiated, the 
unzipping process transforms the stressed crease into two 
overlapping graphene edges in an exothermic reaction. 

Having studied the microscopic cleavage process in 
the warped graphene model, we next turn to hydrogen- 
induced disintegration of specific fullerene and nanotube 
structures. As we discussed above, we expect hydrogen- 
induced disintegration to proceed spontaneously only in 
very large fullerenes, which are typically mult i- wall struc- 
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FIG. 3: Hydrogen- assisted disintegration of (a) the C240 
fullerene and (b) the (6,6) carbon nanotube. Structural snap 
shots along the optimum reaction pathway are depicted in the 
left panels. The structural transformations of the nanotube 
in (b) are shown in end-on view. Carbon atoms are shown 
as grey, hydrogen atoms as larger black spheres. The right 
panels show the total energy change in the system, with the 
data points corresponding to the structures in the left panels. 



tures. In the smaller C240 fullerene, we find this reaction 
to be energetically activated. The energetics of the C240 
disintegration process, induced by hydrogen adsorption 
along the icosahedral edges, is shown in Fig. [3Ka) along 
with snap-shots of intermediate structures. Unlike in the 
model system addressed in Fig. [21 the hydrogen coverage 
has been gradually increased as the cleavage progressed, 
as an alternative scenario to initiating unzipping. Both 
in the warped graphene model and in the C240 fullerene, 
we find the disintegration process to be exothermic, fol- 
lowing an initial energy investment of ~2 eV. 

Whereas large fullerenes are more susceptible to 
hydrogen-induced fracture due to the stress build-up 
in the hydrogenated edges, we find nanotubes that are 
narrow to be more susceptible to fracture due to the 
stress caused by their high uniform curvature. Both in 
fullerenes and nanotubes with multiple walls, hydrogen 
will first adsorb on the outermost wall, causing local de- 
cohesion leading to bond breaking and eventual disinte- 
gration of the outermost wall. Continuing exposure to 
hydrogen will cause a layer-by-layer peeling that will ap- 
pear as exfoliation. 

As suggested by our results in Fig.[T](b), hydrogen ad- 
sorption concentrates the deformation along lines par- 
allel to the tube axis, reducing the flexural stress in 
the hydrogen- free region. The right panels of Fig. [IJb), 
depicting hydrogen covered nanotubes, illustrate how 
hydrogen-induced decohesion along adsorption lines may 
cause large structural changes. The detailed cleavage 
process of an infinitely long (6, 6) carbon nanotube, cov- 
ered by a complete chain of hydrogens adsorbed exohe- 
drally and a neighboring chain adsorbed endohedrally, is 
described in Fig. [3](b). The sequence of snap shots, pre- 
sented in an end-on view, depicts intermediate structures, 
including a scroll [23], encountered during a conjugate 



gradient optimization from an unrelaxed initial structure. 
This process is only spontaneous if all hydrogens adsorb 
in the optimum way to form completely filled lines of ad- 
sorbed hydrogen, and is activated otherwise. The final 
structure with a horseshoe-like cross-section, stabilized 
by the attraction of the hydrogenated zigzag edges, may 
furthermore undergo an activated conversion to a planar 
graphene strip with hydrogen-terminated edges. 

In conclusion, we used ab initio density functional 
calculations to study hydrogen-induced disintegration of 
single- and mult i- wall carbon fullerenes and nanotubes. 
Our results indicate that hydrogen atoms preferentially 
chemisorb along lines in sp^ bonded carbon nanostruc- 
tures, locally weakening the carbon bonds and releas- 
ing stress. Contrary to common belief that chemisorbed 
hydrogen causes no harm to carbon nanostructures, we 
show that irreversible structural changes may be caused 
by chemisorbed hydrogen that can initiate local cleavage 
leading to the disintegration of the outermost wall. In 
experimental observations of multi-wall structures, this 
could be interpreted as hydrogen-induced exfoliation. 
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mance Computing Center. 
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